Experimental evidence for phase synchronization transitions in human 

cardio-respiratory system 
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Transitions in the dynamics of complex systems can be characterized by changes in the synchro- 
nization behavior of their components. Taking the human cardio-respiratory system as an example 
and using an automated procedure for screening the synchrograms of 112 healthy subjects we study 
the frequency and the distribution of synchronization episodes under different physiological condi- 
tions that occur during sleep. We find that phase synchronization between heartbeat and breathing 
is significantly enhanced during non-rapid-eye-movement (non-REM) sleep (deep sleep and light 
sleep) and reduced during REM sleep. Our results suggest that the synchronization is mainly due 
to a weak influence of the breathing oscillator upon the heartbeat oscillator, which is disturbed 
in the presence of long-term correlated noise, superimposed by the activity of higher brain regions 
during REM sleep. 

PACS numbers: 05.45.Xt, 87.19.Hh, 87.19.Uv 



Periodic events are ubiquitous in many natural sys- 
tems [l[ . If two oscillatory processes are weakly coupled, 
they can become phase synchronized. Transitions in the 
synchronization behavior have been shown to be impor- 
tant characteristics of coupled oscillatory model systems 
Q ■ It has also been found that noise, when applied iden- 
tically to different nonlinear oscillators, can induce, en- 
hance, or destroy synchronization among them 
However, phase synchronization is difficult to study in ex- 
perimental data which are very often inherently nonsta- 
tionary and thus contain only quasiperiodic oscillations. 
Among the few recent experimental studies are coupled 
electrochemical oscillators laser systems and cli- 
mate variables Q. In physiology, the study of phase 
synchronization focusses on cardio-respiratory data (see 
below) and encephalographic data [7(. Here, in order 
to obtain reliable experimental evidences of transitions 
in phase synchronization behavior, we consider cardio- 
respiratory synchronization in humans during sleep, be- 
cause homogeneous long-term data for well defined con- 
ditions of a complex system is available in this particular 
example. 

First approaches for the study of cardio-respiratory 
synchronization have been undertaken by the analysis 
of the relative position of inspiration within the corre- 
sponding cardiac cycle [8| . More recently, phase synchro- 
nization between heartbeat and breathing has been stud- 
ied durin g w akefulness using the synchrogram method 
0> 0) II > 12] ■ While long synchronization episodes were 
observed in athletes and heart transplant patients (sev- 
eral hundreds of seconds) [13, EH , shorter episodes were 
detected in normal subjects (typical duration less than 
hundred seconds) [ll|, [lj, [HI ■ For two recent models of 
cardio-respiratory synchronization, see fl4. 15|. 



In this Letter we use the concept of phase synchroniza- 
tion to develop an automated synchrogram based proce- 
dure and study interactions between cardiac and respira- 
tory oscillations under different well-defined physiological 
conditions. We focus on the sleep stages, where external 
stimuli are absent. It is well known that healthy sleep 
consists of cycles of roughly 1-2 hours duration. Each 
cycle is characterized by a sequence starting usually with 
light sleep, followed by deep sleep and REM sleep (rapid 
eye movement) [16]. We find the intriguing result that 
during REM sleep cardio-respiratory synchronization is 
suppressed by approximately a factor of 3 compared with 
wakefulness. On the other hand, during non-REM sleep, 
it is enhanced by a factor of 2.4, again compared with 
wakefulness. In addition, we find that these significant 
differences between synchronization in REM and non- 
REM sleep are very stable and occur in the same way 
for males and females, independent of age and indepen- 
dent of the body mass index (BMI) . Hence it seems likely 
that - similar to the long-term correlations [l7| occurring 
in both heartbeat [r| and breathing [r| during REM 
sleep but not during non-REM sleep - the differences are 
caused by the influence of the activity of higher brain 
regions on both oscillators [20j ]. 

First, we developed an algorithm, which detects 
epochs of synchronization automatically and systemat- 
ically. The algorithm is applied on simultaneous records 
of respiration (from a thermistor placed close to the sub- 
ject's nose) and heartbeat (from electrocardiograms) ob- 
tained for 112 healthy subjects during sleep. The data 
was recorded in the EUproject SIESTA in several Euro- 
pean sleep laboratories [21| , and the average length of the 
records is 7.9 hours with a standard deviation of 25 min- 
utes. Sleep stages have been determined by visual evalu- 
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FIG. 1: Illustration of our automated phase synchronization 
detector: (a) section of a typical synchrogram in raw form, 
(b) the same section as in (a) after applying a moving average 
filter with windows of r = 30 s around each breathing cycle; 
the window averages are shown together with their standard 
deviations (error bars), (c) contiguous phase points of (b) 
where the standard deviation is below a threshold (see text) 
are detected as synchronization episodes if their durations 
exceed T seconds (here T = 30 s). 



ation of electrophysiological recordings of brain activity 
We can thus assign the synchronization episodes 
to the specific sleep stages. In addition, we constructed 
surrogate data by random combination of heartbeat and 
breathing signals from different subjects We found 
that the total duration of the detected synchronization 
episodes in real data is increased by a factor 2.3 as com- 
pared with the surrogate data, suggesting that most of 
the detected episodes are real 



22] 



Our algorithm for the detection of phase synchroniza- 
tion episodes is based on the study of cardio-respiratory 
synchrograms [sLGiJ 11 , 12 1 . For each record, the times tk 
of heartbeats are mapped on the continuous cumulative 
phase $ r (t) of the respiratory signal, which we obtain by 
unfolding the phase of the Hilbert transform of the nor- 
malized thermistor recording. Figure 1(a) shows a repre- 
sentative synchrogram, where <fr r (tk) mod 27r is plotted 
versus tk- In case of n:l synchronization (i. e., if n heart- 
beats fit to one breathing cycle) one observes n parallel 
horizontal lines in the synchrogram [n — 3 in Fig. 1(a)]. 
In general, to find different ratios n:m of phase synchro- 
nization, we plot \& m (tfe) = $ r (tfc) mod 27rm versus tk- 

While most earlier work relies on a visual evaluation 
of the synchrograms 

nana, 

we detect the episodes in 
a fully systematic way. For each synchronization ratio 
n:m we first replace the n phase points *$> m (tk) in each 
m respiratory cycles by the averages ^ m {tk) calculated 
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FIG. 2: Medians, upper and lower quartiles (bars) and means 
(filled symbols) for the detected synchronization rates (a) ver- 
sus T for all original data (dotted bars and circles) and sur- 
rogate data (striped bars and triangles), (left of dotted line 
5 = 6, right of dotted line 8 = 5); (b,c) versus age, body mass 
indices (BMI) and gender, for wakefulness (dotted bars and 
circles), REM sleep (blank bars and triangles) and non-REM 
sleep (striped bars and diamonds) for T = 30 s. Note the 
similar synchronization behavior in all groups. 



over the corresponding points in the time windows from 
tk — t/2 to tk + t/2 [Fig. 1(b)]. In the second step, the 
algorithm deletes all phase points \& m (tfc) where the mean 
standard deviation of the n points in each m breathing 
cycles, (cr) n , is larger than 2rmr/n6. In the third step, 
only the phase points \& m (ifc) in uninterrupted sequences 
of durations exceeding T seconds are kept [Fig. 1(c)]. 

Figure 2(a) shows a comparison of the detected syn- 
chronization rates in real data and in surrogate data 
for several values of the parameter T. The ratio of the 
mean percentage of synchronization in real data over the 
mean percentage in surrogate data increases from 1.6 for 
T = 20 s to 3.4 for T = 40 s. There is no generic limit 
for T. However, we choose T = 30 s in order to keep 
the number of arbitrary detection of synchronization at 
an acceptable low level (ss 1% in the surrogate data) 
while at the same time still detecting synchronization 
episodes during all kinds of sleep stages. In addition, for 
T = t = 30 s we have only one effective time scale pa- 
rameter, which is identical with the standardized time 
frame used for the detection of sleep stages [lj| • We note 
that changing the parameter S has a similar effect on the 
results as changing T, and we chose 6 = 5 based on simi- 
lar considerations as for T. Our results do hardly depend 
on the duration r of the initial running average. 

When studying the percentages of synchronized 
episodes in real data separately during wakefulness, REM 
sleep, and non-REM sleep we obtain a highly significant 
difference in the frequency of cardio-respiratory synchro- 
nization between the two major physiological states dur- 
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ing sleep. We find 3.8% synchronization in non-REM 
sleep compared with just 0.6% in REM sleep - a differ- 
ence by a factor of 6.3. Wakefulness is clearly interme- 
diate, since we find 1.6% for it. Similar differences are 
observed for other values of T and 6. Since our data 
base contains records of 112 subjects, we can also study 
synchronization separately for several age groups, several 
groups with different body mass index (BMI), and men 
and women. Figures 2(b,c) show that the results, for 
wakefulness, REM sleep, and non-REM sleep are prac- 
tically the same for both genders, all BMI groups, and 
all age groups, although both, heart rate and breathing 
rate are known to depend on BMI and age. These results 
prove that our finding of significant differences between 
the cardio-respiratory synchronization in REM and non- 
REM sleep is very stable. 

Similar stable differences between REM and non-REM 
sleep were found in the correlation properties of both 
heartbeat [3] and breathing [H| fluctuations, but hardly 
in the magnitude of these fluctuations |20j. The differ- 
ences were attributed to the influence of the central ner- 
vous system with its sleep stage regulation in higher brain 
regions on the autonomous nervous system. The similar- 
ity leads us to suggest that the diminished synchroniza- 
tion during REM sleep is also caused by influences of 
the central nervous system. As long as the heartbeat 
oscillator and the breathing oscillator (as parts of the 
autonomous nervous system) are only affected by uncor- 
rected noise from higher brain regions, they run like two 
weakly coupled oscillators - and they clearly show syn- 
chronization as expected, possibly enhanced by the noise 
0| . However, if the higher brain regions are more ac- 
tive and impose long-term correlated noise on the two 
oscillators, as is the case during REM sleep, the noise 
disturbs the emergence of synchronized patterns, lead- 
ing to a drastic reduction of synchronization episodes. 
Hence we suggest from the experimental data that cor- 
related noise is rather suppressing synchronization while 
uncorrelated noise might increase it. 

Our interpretation is consistent with the result that 
cardio-respiratory synchronization is enhanced in heart 
transplanted patients, where correlated signals from the 
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brain can hardly affect the heartbeat oscillator 
Hence, it supports that any relation of the synchroniza- 
tion patterns with cardiac impairments can only be an 
indirect one as reported recently [2~3| . Diminished long- 
term correlated regulation activity might explain the in- 
crease of synchronization in well-trained athletes [loj ] . 
where fluctuations of heartbeat and breathing might be 
avoided to optimize the cardiovascular system for opti- 
mal performance. 

In order to gain insight into the mechanism of cardio- 
respiratory phase synchronization we have studied, again 
for all 112 subjects, the distribution of the synchroniza- 
tion ratios n:m, where n cardiac cycles are synchronized 
with m breathing cycles. Figure 3 shows the normal- 
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FIG. 3: Normalized histograms of all detected synchroniza- 
tion ratios as well as distributions of the frequency ratios be- 
tween heartbeat and breathing during (a) wakefulness, (b) 
REM sleep, (c) light sleep, and (d) deep sleep. Note that 
4:1 synchronization is remarkably increased in light sleep and 
deep sleep, whereas in REM sleep 3:1 synchronization dom- 
inates. In (e) the quotient of the synchronization ratio his- 
togram with the distribution of the frequency ratios is shown 
for the original data (crosses) as well as for the surrogate data 
(circles). 



ized histograms of the synchronization ratios during (a) 
wakefulness, (b) REM sleep, (c) light sleep, and (d) deep 
sleep. The underlying continuous curves in Fig. 3(a- 
d) show the distributions of the frequency ratios between 
heartbeat and breathing independent of synchronization. 

To clarify the efficiency of the synchronization mech- 
anism, we show the quotient of the synchronization rate 
histogram with the distribution of the frequency ratios 
(crosses) in Fig. 3(e). When comparing with the corre- 
sponding curve for surrogate data (circles) it becomes ob- 
vious that n:l synchronization is preferred. This involves 
not just the common ratios of 3:1, 4:1, and 5:1, but also 
clearly 6:1 synchronization. In particular, we do not find 
any indication of a suppression of 4:1 synchronization as 
has been reported recently based on modelling [l4] . This 
suggests that the feedback from baroreceptors to the res- 
piratory centers introduced in the model is probably not 
very important in the healthy subjects we studied. The 
synchronization ratios n:2 and n:3 are weakly efficient 
if n is low, but not efficient at all for large values of n. 
Hence, Fig. 3(e) proves that cardio-respiratory synchro- 
nization is nearly limited to frequency ratios n:m with 
very small m, but quite independent of the number of 
heartbeats n. 

^From this behavior we suggest that the physiological 
synchronization mechanism is mainly based on an inter- 
action of the respiratory cycle upon heartbeat and not 
vice versa. This assumption can explain the weak (or 
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absent) efficiency of n:m synchronization with m > 1. 
If the heartbeat oscillator gets a synchronizing kick at 
a particular phase of each respiratory cycle and m > 1, 
only half or even less of the kicks coincide with a heart- 
beat and thus n:m synchronization cannot be effective. 
The assumption is consistent with the result obtained 
when studying the direction of synchronization in chil- 
dren and adults Q and in a recent model [H| . It is also 
coherent with the result that synchronization is enhanced 
under paced respiration (l3| . 

In conclusion, we have studied cardio-respiratory phase 
synchronization during different well-defined physiolog- 
ical stages in sleep for a large data base of healthy 
subjects. We observed clearly reduced synchronization 
during REM sleep and enhanced synchronization during 
non-REM sleep. The result is stable for all studied sub- 
groups of subjects; it is neither affected by gender, nor by 
age, nor by BMI. Since REM and non-REM sleep differ 
mainly in the type of activity of higher brain centers, it 
seems probable that the differences in cardio-respiratory 
synchronization are caused by the more and less long- 
term correlated regulation actions of the brain during 
REM and non-REM, respectively. Heart rate and breath- 
ing rhythm generators behave like two weakly coupled os- 
cillators, where the coupling direction is from breathing 
to heartbeat. They become synchronized if uncorrelated 
noise is imposed from the brain while long-term corre- 
lated noise disturbs the emergence of the synchronized 
patterns. Hence, the experimental data suggests that 
correlated noise is suppressing synchronization while un- 
correlated noise might increase it. 
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